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Based on recent results obtained by the Pierre Auger Observatory (PAO), it has been hypothesized 
that Centaurus A (Cen A) is a source of ultra-high-energy cosmic rays (UHECRs) and associated 
neutrinos. We point out that the diffuse neutrino flux may be used to constrain the source model 
if one assumes that the ratio between the UHECR and neutrino fluxes outputted by Cen A is 
representative for other sources. Under this assumption we investigate the relation between the 
neutrino flux from Cen A and the diffuse neutrino flux. Assuming furthermore that Cen A is the 
source of two UHECR events observed by PAO, we estimate the all-sky diffuse neutrino flux to be 
~200 — 5000 times larger than the neutrino flux from Cen A. As a result, the diffuse neutrino fluxes 
associated with some of the recently proposed models of UHECR-related neutrino production in 
Cen A are above existing limits. Regardless of the underlying source model, our results indicate 
that the detection of neutrinos from Cen A without the accompanying diffuse flux would mean that 
^ ' Cen A is an exceptionally efficient neutrino source. 
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P-C I. INTRODUCTION 

i 

The Pierre Auger Observatory (PAO) has recently reported new results [H, on the arrival directions of the highest- 
energy cosmic rays (CRs). The data show strong evidence for anisotropy of these CRs, which suggests that at least 

1 i ' some sources of ultra- high-energy cosmic rays (UHECRs) are relatively close. The cosmic-ray anisotropy has been 

confirmed by other studies using different statistical methods and source catalogs 0, 0, S] ■ The correlation between 
the arrival directions of these CRs and the positions of known active galactic nuclei (AGNs) has lead the PAO to 
suggest that nearby AGNs, or astrophysical objects with a similar spatial distribution, are the sources of UHECRs 
However, the observed deficit of UHECRs from the nearby Virgo cluster appears to be incompatible with such a 
source distribution [fj]. Therefore the origin of UHECRs remains unclear at present. 
■ Independent of a possible general connection between UHECRs and AGNs, the PAO data raise the possibility that 
Centaurus A (Cen A) is a source of UHECRs. The PAO analysis [l|, associates two events with Cen A, but it has 
been pointed out that at least four events can be associated with Cen A if one takes account of its morphology 0]. 
' At a distance of ^3.5 Mpc, Cen A (NGC 5128) is the nearest active galaxy (see Ref. [§[ for a review). It is classified 
as a Fanaroff-Riley type I radio galaxy, possibly harboringa misdirected BL Lac nucleus. The galaxy is believed to 
be powered by accretion on a ^10 8 Mq central black hole [9(. It has a very compact nucleus, a pronounced northern 
jet, a dimmer southern jet, and giant radio lobes extending out to ~250 kpc. We likely observe the jets from the side, 
under a viewing angle of 50° — 80° [l(| (see, however, Ref. [Hj]). Due to its proximity, it has been suggested long ago 
that this galaxy may be the source of UHECR events observed at Earth [12, El, E3 • 

A general prediction of models of UHECR acceleration is the production of neutrinos through the interaction of 
accelerated protons (or nuclei) with the ambient photon field or with target protons in the source. Hence, for a given 
acceleration mechanism the neutrino and UHECR fluxes are related. Detection of neutrinos from Cen A or limits on 
their flux may therefore translate into constraints on the underlying acceleration models. 

Neutrino production in AGNs has been studied by many authors, see e.g. Refs. [H El, El El, El, US HI HI HI, 
24], HI]. More recently, the authors of Refs. [H,n3] have presented estimates on the neutrino flux from Cen A under 
the assumption that two out of the 27 UHECR events in the PAO analysis can be attributed to this galaxy. These 
authors have however not considered the diffuse flux due to all (unresolved) neutrino sources within their models. 
Assuming that the environment in Cen A is somehow representative for all sources of UHECRs and accompanying 
neutrinos, the diffuse neutrino flux may also constrain the source model. In this work we investigate the connection 
between the UHECR-related neutrino flux from Cen A and the associated diffuse neutrino flux. 
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In relating the diffuse flux to that from Cen A, we assume (as a working hypothesis) that Cen A is a 'typical' 
source of UHECRs and neutrinos, i.e. we assume universal UHECR and neutrino injection spectra with a fixed 
relative strength. We do not make any assumptions on the intrinsic luminosity or on the distance of the UHECR 
sources. Within the assumption of typicality, the diffuse neutrino flux can be estimated by upscaling the neutrino 
flux from Cen A using CR data. The scaling factor is the product of a trivial factor standing for the fraction of 
observed UHECR events that is attributed to Cen A, and a non-trivial factor that accounts for the difference in CR 
and neutrino mean free paths: as the UHECR flux from far-away sources is strongly attenuated by interactions with 
the cosmic microwave background, the ratio of the diffuse neutrino flux to the Cen-A neutrino flux will be much larger 
than the observed ratio of the diffuse UHECR flux to the Cen-A UHECR flux. In this study we estimate this scaling 
factor without relying on any specific source model. For definiteness we assume that the UHECRs are protons, the 
composition of UHECRs still being under debate [28| . 

We find that the all-sky diffuse neutrino flux is expected to be ^200 — 5000 times larger than the neutrino flux 
from Cen A, depending most strongly on the assumed source evolution. As a consequence, diffuse neutrino fluxes 
associated with some of the recently pro posed models of UHECR-related neutrino production in Cen A overshoot 
existing bounds by the AMANDA-II [29l . l30l ] and PAO [3l| experiments. Regardless of the underlying production 
model, our results indicate that the detection of neutrinos from Cen A without the accompanying diffuse flux would 
imply that Cen A is an exceptionally strong neutrino source. 

This paper is organized as follows. In section [II] we discuss attenuation of the UHECR proton and neutrino fluxes 
and estimate the ratio of the diffuse neutrino flux to the Cen-A neutrino flux. In section UHl we apply these results 
to models that were recently proposed in Refs. [H, [13] • We summarize and discuss our work in section [TVl The 
appendices contain additional information on the computer code used to calculate attenuation of the cosmic proton 
flux, and an estimate of the neutrino effective area of the IceCube experiment. 



II. RELATING THE DIFFUSE PROTON AND NEUTRINO FLUXES 



The flux of high-energy protons from a cosmic source is attenuated by redshift and by interactions with CMB 
photons. As a result, the observed flux of UHECR protons is significantly smaller than the flux that is injected by 
all sources. Neutrinos, on the other hand, only suffer redshift energy losses, which is of far lesser importance. This 
difference boosts the diffuse neutrino flux reaching Earth compared to the diffuse flux of UHECRs, an effect that 
should be taken into account when normalizing the diffuse neutrino flux to the observed UHECR flux. Under the 
assumption of a universal relation between the outputted neutrino and proton fluxes, this effect may be parameterized 
by a parameter H that we introduce in this section. 



A. The neutrino boost factor 



The observed differential flux <f> from a single source at proper distance D (redshift z) can be expressed as 

j°(E ) dE 



<f>{E) 



47rD 2 (l + z) dE 



(1) 



where E is the observed energy, Eq = Eq (E, z) is the energy at the source, and j° denotes the differential spectrum 
at the source. Integrating over a cosmological distribution of sources, the diffuse flux is equal to (see, e.g., Ref. [33 |) 
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where no is the present source density and e(z) parameterizes source evolution (no evolution corresponds to e{z) = 1) 
In this expression we assume that all sources are identical. Within the ACDM concordance model that we adopt, 



1 



h (i + 2)va»(i + 2) 3 + ^A 



(3) 



where Hq — 73 km s" 1 Mpc -1 denotes the present Hubble constant, Sl m = 0.24 is the present matter density 
parameter, and J7a = 0.76 is the present vacuum energy density parameter [33|. 

Since we assume a universal neutrino injection spectrum in this work, we approximate the observed diffuse neutrino 
flux by <f>f s (E) oc £(E). (In the case of spectral breaks this neglects smearing due to different source redshifts. We 
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will come back to this issue in the next section.) We normalize the diffuse neutrino flux to the integral UHECR flux 
$p above a threshold energy E t h, and define a constant of proportionality H as follows: 

where j® is the differential neutrino spectrum at the source, and J® is the integral UHECR proton spectrum at the 
source (we will use capital symbols to refer to integral spectra and fluxes, and lower-case symbols for differential ones). 
The effect of the different UHECR proton and neutrino mean free path lengths is now contained in the scaling factor 
H, which we will refer to as the neutrino boost factor. It can be determined for any neutrino and any proton injection 
spectrum from eq. together with a formula for Eq(E, z). Note that H depends on the threshold energy E t h but 
not on the neutrino energy E. 

We have computed the parameter H numerically for power-law proton and neutrino injection spectra, jp(E) cx E~ Pp 
and fl (E) cx E~ Vv , respectively. This is done with the help of a computer code that solves eq. ^ for protons and 
neutrinos and then determines H from eq. Proton energy losses are taken into account in the continuous loss 
approximation using expressions for the energy-loss time that are given in appendix [X] In this process we integrate 
over a cosmological distribution of sources up to redshift z = 5. Because we do not know the redshift evolution of 
UHECR sources, we consider as limiting cases both no evolution and strong evolution tracing the AGN luminosity 



density evolution given in Ref. 34[, i.e. e(z) cx (1 + z) 3 4 up to z = 1.9, a constant e up to z = 3, and negative 
evolution e(z) cx (z — 3)~ 033 beyond z = 3. We have not explicitly considered milder source evolution models, but 
the results for H will be numerically intermediate between the two cases considered. 

In figure [T] we show the neutrino boost factor H as a function of proton power-law index p p and for different 
neutrino power-law indices p u . In producing this figure, we have taken the detector threshold energy equal to Eth 
= 57 EeV (the energy threshold used in the PAO analysis [3, 0), and we have assumed that the maximum proton 
energy is much larger than this. As can be seen in the figure, the boost factor increases mildly as the neutrino 
injection spectrum becomes harder compared to the proton injection spectrum, which can be understood from the 
scaling of energy spectra with redshift. Including source evolution is a more dramatic effect, increasing H by an order 
of magnitude. This of course reflects the fact that, with strong source evolution, the fraction of sources that can be 
seen in high-energy neutrinos but not in UHECRs increases substantially. 

The systematic uncertainty of ^20% in energy determination by PAO f35j introduces some uncertainty in our 
results on the neutrino boost factor H . If the actual threshold energy is lower than 57 EeV, the effect of UHECR flux 
attenuation is less severe and hence the neutrino boost factor H is smaller than the results presented in figure [1] We 
have found that, for the parameters used in figure [I] this effect is well approximated by the simple formula 

H(E th ) = lO^T- 1 H(E Q ), (5) 

where E t h now denotes the actual threshold energy and Eo = 57 EeV. As this equation shows, the neutrino boost 
factor H becomes smaller (larger) by a factor ~1.6 in case the energy is systematically overestimated (underestimated) 
by 20%. 



B. Scaling the neutrino flux from Cen A 



We now specialize the discussion to the recent PAO data. The number of events with energy above Eth from a 
point source at declination S s can be expressed as A pt = & pt (E t h)TA(5 s ), where <& pt (£ , t h) is the integral flux above 
E t h, T is the observation time, and A(S S ) denotes the experiment's effective area for a source at declination 8 S . The 
number of events due to the diffuse flux is then N dlS = <j> dlfT E!, where <E> dlff is the integral diffuse flux per sterad, and 
S = T J A(5 s )dil denotes the exposure. With these expressions, we estimate the diffuse UHECR flux $^ and the 
UHECR flux from Cen A ^ cnA above E th = 57 EeV to be: 

$ diff (E th ) = — ~j VconA = 9 x 10- 21 cm- 2 s-hr- 1 ; (6) 

^ A (Eth) = N ^SMSs)dsi = g x 10 - 2Ws -! (7) 

s A(d s ) 

where S = 9000 km 2 yr s is the total Auger exposure, S s = —43° is the declination of Cen A, iV tot = 27 is the 
total number of observed UHECR events, and iVcenA is the number of events from Cen A. Following the PAO 
analysis [l], HJ we attribute Ac C nA = 2 events to Cen A (note however that the actual number may be larger 0]). 
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FIG. 1: Neutrino boost factor H as a function of proton power-law index p p for five neutrino power-law indices p v and for 
£th = 57 EeV. From top to bottom the neutrino spectrum becomes progressively softer compared to the proton spectrum: 

Pv = Pp — 1.0 (dashed line), p v — p p — 0.5 (dash-dash-dotted), p v = p v (solid), p„ = p v + 0.5 (dot-dot-dashed), and 
Pv — Pp + 1.0 (dotted). In the left panel we assumed no source evolution; in the right panel we assumed that the sources 

follow AGN evolution. 



In eq. we used the relation A(S S ) oc u)(S s ), where u> is the relative PAO exposure given in Ref. [3(|, to estimate 
A(S S )/ f A(5 s )dn = 0.15 sr" 1 . 

Using the definition (j4| of the neutrino boost factor H, we now express the ratio of the neutrino flux from Cen A 
to the associated diffuse flux as follows: 

<j>**(E) H(E th )$£ s (E th ) H(E th ) (N tot - N CcnA )A(6 s )_ = x g^^-i ^ (g) 



0Cc„A ( £) $CenA (i ; th) Nc ea A J A(S s )dQ 

where we have used the fact that protons from Cen A reach Earth virtually without energy loss. This equation is the 
main result of this paper. Together with the results presented in figure ^ it allows to estimate the diffuse neutrino flux 
from a model neutrino flux for Cen A, under the assumption that the physical environment in Cen A is representative 
for all UHECR and neutrino sources. Within the parameter range shown in figure [TJ we thus find that the all-sky 
diffuse neutrino flux is ^200 — 500 (800 — 5000) times larger than the neutrino flux from Cen A in the case of no 
(strong) source evolution. 

We now compare the expected event rate in a neutrino detector for neutrinos from Cen A to the event rate for the 
diffuse neutrino flux. In this work we focus on neutrino detection with IceCube for definiteness. Neutrino detection 
with IceCube is discussed in appendix [5] Neutrinos from Cen A are downgoing for IceCube. Hence, for a fair 
comparison between Cen A and the diffuse flux, we consider both upgoing and downgoing diffuse neutrinos although 
detection of the latter is complicated by the atmospheric muon background. From eqs. ®, (|B6[) . and (|B8|) we find 
that the ratio of events associated with the diffuse downgoing neutrino flux to events associated with Cen A is 

jydiff , dn 

^a- 11 ^), (9) 

where we took the IceCube field-of-view equal to fij = 5.7 (up to 5° from the horizon). In this expression we 
approximated the effective area for downgoing neutrinos as angle-independent. The ratio of events associated with 
the diffuse upgoing neutrino flux to events associated with Cen A is 

Ardiff, up 

llx#(£th), (10) 



where (cf. eqs. (|B6j and (|B9| ) 



jVCcnA 



_ JdEME)AZ s (E) 
X - JdEME)A% s (E) [ > 



is a factor of order unity that accounts for the angular dependence of the detector. Here A^ n cS (A^ p cS ) denotes 
the (average) effective area for downgoing (upgoing) diffuse neutrinos. Using estimates for the effective areas of the 
IceCube detector given in eqs. (|B7[) and (|B10|) . respectively, we have verified that 1 < % < 2 for neutrino test spectra 
<f> v oc E-p, where 1 < p < 3 and 10 3 GeV < E v < 10 s GeV. 
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III. SOURCE MODEL EXAMPLE 



Following the hypothesis that 2 of the 27 UHECR events detected by PAO are from Cen A, different models for 
ultra-high-energy (UHE) proton and neutrino production in Cen A were proposed in Refs. pq . [2?| • For reasons of 
space we focus on the model adopted in Ref. [26J, hereafter referred to as the CH model, which is based on earlier 
work in Ref. [l8| . For this model, we estimate the associated diffuse neutrino flux and compare its detection prospects 
to those of the neutrinos from Cen A using the general results obtained in the previous section. In the last part of 
this section we comment on the detection prospects of the models that were recently proposed in Ref. j2?| ■ 



A. The model 

An attractive feature of the model adopted in Ref. [26| is that the neutrino flux at high energies is harder than the 
flux of UHECR protons. To achieve this, the model requires that a population of high-energy seed protons (accelerated 
through e.g. shock acceleration) is confined to a region close to the source. These protons create neutrons and neutrinos 
in photopion interactions with the ambient photon field. The neutrons escape from the source, decay, and give rise 
to UHECR protons. In this process the neutrons lose energy in interactions with the photon field before decay, thus 
softening the spectrum of UHE protons. Neutrinos, on the other hand, trace the neutron energy spectrum upon 
production (i.e. without energy loss). 

The model predicts two spectral breaks in the cosmic-ray injection spectrum at the energies where the optical 
depths for proton and neutron photopion production become unity (see Ref. fl8j for a thorough discussion). Since 
these energies are generally similar, one may assume a single break energy E^ t . Below this break energy both the 
cosmic-ray injection spectrum and the neutrino spectrum are harder than the seed proton spectrum by one power of 
the energy (assuming photopion production on a n 7 (e 7 ) oc e~ 2 photon field). Above the break energy the cosmic-ray 
injection spectrum is softer than the seed proton spectrum by one power of the energy (assuming the same photon 
field) , while the neutrino spectrum follows the initial seed proton spectrum. Identifying the high-energy part of the 
injected cosmic-ray spectrum with the observed UHECR flux <p p , we can express the all- flavor neutrino flux </>„ x as: 

&, . ( E E 2 \ / E 



^ (s) = £^ min l^^'^^J^ \^r) ' (12) 

Zn'lvn \'/im^br 'hn^bi / \'\vn J 

where £j is the energy fraction of the proton that is transferred to species i (neutron or neutrino) in photopion 
interactions, and r\ vn is the ratio of the average neutrino energy to the average neutron energy. The values of these 
quantities depend on the spectral distribution of the photon field and can be estimated numerically. The authors of 
Ref. [2j| take £j//£ n — 0-2 and r\ vn = 0.04 for interactions in the nucleus of Cen A. The break energy, being due to 
a change in photopion production efficiency, is determined by the ambient photon distribution. Although its value 
cannot be directly inferred from observations, it may be estimated from the observed gamma-ray spectrum because 
interactions between gamma rays and the ambient photons also give rise to a spectral break in the gamma-ray flux. 
In this way the break energy in the UHECR spectrum can be estimated as 2?br ~ 3 x 10 s -E 7 ,br, where i£y,br denotes 
the gamma-ray break energy [l8|, [26j|. The authors of Ref. [26| conservatively take -E 7 .br — 200 MeV for Cen A, so 
that £ b r ^ 10 8 GeV. 

In this work we only consider detection of muon neutrinos for definiteness, and hence we should take account of the 
neutrino flavor ratios. In the CH model, the source is optically thin and the effect of meson synchrotron energy loss 
is neglected. Under these conditions, the neutrino flux from proton-photon interactions is dominated by neutrinos 
from pion decay at all energies (37j . This implies that the neutrino flavor ratio (electron : muon : tau) at the source 
is approximately 1 : 2 : 0, as is also indicated in fig. 1 of Ref. [26|]. In this case neutrino oscillations lead to a flavor 
ratio close to 1 : 1 : 1 at Earth [38[ . Hence the expected muon- neutrino flux from Cen A, in the CH model, is: 



3 £n Eth \ Eth J \ E„ t br / \ E v . b 

3 x 10_n f 77^77 ) min (l, -J^—\ GeV -1 cm -2 s -1 , (13) 
VlCev/ V E v , br/ 



where E v ^ t = r\ vn E\,x — 4 x 10 6 GeV. 
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FIG. 2: Diffuse muon-neutrino flux obtained by scaling (thick solid line) and by numerical computation (solid), together 
with the diffuse UHECR flux obtained numerically (dotted) and the atmospheric muon-neutrino background (dashed) as a 
function of neutrino energy. Also shown are existing upper limits from AMANDA-II and PAO, the projected upper limit for 
IceCube (3 years), and UHECR data from PAO (diamonds). Left panel: no source evolution; right panel: strong source 

evolution. 



B. Diffuse neutrino flux 



The neutrino boost factor for a proton power-law spectrum with index p — 2.7 and a neutrino power-law spectrum 
with index p — 1.7, with and without source evolution, is H = 19 and H = 159, respectively (see fig. [TJ. Scaling 
the model neutrino flux from Cen A from eq. (fT3"|) , we find the following estimates for the diffuse muon-neutrino flux 
within the CH model: 

*»<*> * 10- (ji^Y" min (l,^-) GeV-Ws- St - (14, 

in the case of no source evolution, and a factor of 10 larger in the case of strong evolution. These fluxes are shown in 
figure [2J Also shown are numerical results (obtained by the method described in section [TTJ) for the diffuse UHECR 
and neutrino fluxes after propagation, the preliminary AMANDA-II UHE limit [3(|, the PAO limit [3l|, the projected 
limit for IceCube [3^| after three years of data-taking, PAO data on the UHECR proton flux [4(| , and the atmospheric 
neutrino background. The AMANDA-II, PAO, and IceCube detection limits are the 90% confidence level upper limits 
for a <j) v oc E~ 2 diffuse muon-neutrino flux, where we have assumed a 1 : 1 : 1 flavor ratio at Earth. The atmospheric 
muon-neutrino background flux is parameterized as follows: 

f 8-4xl0- 2 (^/l GcV)- 2 74 r y-l rTT1 -2 _-l „--! ( F , -, n 5.3 CpV \ 

(h h s(E) = l — 1+0.002(^/1 Gov) — (jeV cm s sr IA> < w UeV J H5) 

VvKv) \ 5.7 x 10- 3 (^/lGeV)" 3 ' 01 GeV" 1 cm" 2 s" 1 sr" 1 (E v > 10 5 3 GeV) ' 

where the high-energy contribution is due to prompt charm decay. This parameterization is, by construction, close 
to the maximum background indicated in figure 11 of Ref. 39]. 

We observe from figure [2] that the scaled neutrino flux is an excellent approximation to the numerical results in the 
high-energy regime, where the energy spectrum follows a single power-law. The spectral break in the scaled neutrino 
flux is however much sharper than in the numerical results, where it is smoothened due to an averaging over redshift. 
As a result, the approximation obtained by scaling underestimates the resulting diffuse neutrino flux at lower energies. 
We have verified that this does not strongly affect the expected number of neutrino events in IceCube. (In principle, 
the diffuse flux at energies below and above the break can be estimated independently from the results obtained is 
section UH which would give a better estimate. Given the other uncertainties we are faced with, we will not pursue 
this though.) In producing this figure, we have assumed that the maximum neutrino energy is larger than 10 11 GeV. 
The exact value is not very important for our estimates, as the interaction rate in IceCube is dominated by neutrinos 
of energy below 10 8 GeV. 

We note from the figure that the UHECR flux follows the PAO data above 57 EeV, which is a consistency check 
of our numerical method. The model flux for strong evolution is marginally consistent with the PAO data at lower 
energies, given the large uncertainties in energy calibration. 

The estimated diffuse neutrino flux for the CH model is well above AMANDA-II and PAO limits in the case of 
strong evolution. This implies that either the CH model (when applied to strongly evolving sources) is too optimistic, 
Cen A is intrinsically exceptional, or that significantly more than two of the UHECR events observed by PAO are 
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from Cen A. In the case of no evolution, the diffuse flux is only marginally above these limits. (Note that the limits 
are based on a neutrino flux <\> v oc E~ 2 , so that a more accurate comparison requires further analysis.) The flux is 
however well above the projected IceCube upper limit. At energies larger than ~10 6 GeV, the atmospheric neutrino 
background is strongly suppressed and hence detection should pose no difficulties. 



C. Event rates 



Using estimates for the neutrino effective area presented in appendix [B] (see fig. [3]), we can now estimate the 
neutrino event rate of the diffuse neutrino flux and that from Cen A in the IceCube neutrino detector. Since Cen A is 
in the southern hemisphere, neutrinos from this galaxy are downgoing for IceCube. This makes detection challenging, 
and only possible at very high energies. To eliminate the background we consider the flux of neutrinos with energy 
10 6 GeV < E v < I0 8 GeV (cf. fig. [5]). For a fair comparison between the neutrino interaction rates for Cen A and 
the diffuse flux, we consider both up- and downgoing diffuse neutrinos. We find that x = 1.4 for the CH model, where 
X is the factor that enters the scaling relation eq. (fTU|) . We adopt a field of view £1/ = 5.7, corresponding to cutting 
at 5° below or above the horizon. 

From eq. (|B6[) , we find 0.08 events per year for Cen A with an expected background of 0.004 (6>/5 ) 2 , where 9 is 
the angular resolution. This is in reasonable agreement with results obtained in Ref. [261 ] , where an event rate of 
0.35 per year is found for neutrinos of all flavors (assuming equal detection probabilities). Using the scaling relations 
eqs. © and (jTUJ) , we directly obtain our estimate of 17 (24) neutrino events per year due to downgoing (upgoing) 
diffuse neutrinos per year in the case if no source evolution. Including source evolution, we find 145 (203) downgoing 
(upgoing) events per year. In this energy range the number of background events for the downgoing (upgoing) diffuse 
neutrino flux is roughly 0.4 (1) per year, and hence detection is virtually background-free. 

Our results for the event rate due to the diffuse flux of upgoing neutrinos in the CH model are larger than recent 
estimates in Ref. 41], who estimate ~5 events per year for this model. This result is to be compared with our 
no-evolution estimate of 24 events per year. The authors of Ref. 4l| attribute the diffuse neutrino and UHECR flux 
to Fanaroff- Riley I (FRI) radio galaxies, of which Cen A is an example. They estimate the diffuse neutrino flux by 
adding the contribution of sources up to z = 0.5 using the inferred FRI source density. In contrast, our estimates 
are normalized to the UHECR flux and we consider sources up to z = 5. It is reassuring that these results, which 
are obtained in a different manner, are within a factor few. The fact that our estimates are somewhat larger may be 
attributed to the larger maximum redshift that we have adopted (see also the comment at the end of section III of 
Ref. 111). 



D. Comparison with other models 

In the above we have focused on a model for Cen A that was put forward in Ref. (26[. More recently, the authors 
of Ref. [13] have also discussed UHE proton and neutrino production in Cen A. The authors consider three different 
spectra for the accelerated protons: (i) a straight power-law spectrum with index p = 2.0, (ii) a broken power-law 
spectrum with p — 2.0 before and p — 2.7 after the break energy, and (iii) a straight power-law spectrum with 
index p = 1.2. The first two spectra may be the result of stochastic shock acceleration, the last spectrum of linear 
acceleration in a regular electric field. In case of acceleration near the core neutrinos are produced predominantly in 
the interaction of the accelerated protons with UV photons while low-energy protons provide the dominant target for 
neutrino production if the protons are accelerated in the jet [27] ]. The authors obtain neutrino spectra for the three 
acceleration mechanisms, applied to both the core and the jets, numerically. We do not attempt to reproduce these 
here, but rather estimate the associated diffuse neutrino flux and the event rates from their results. 

For the three acceleration models (i), (ii) and (iii), we find that H = 14, 13, and 14, respectively, in the case of no 
evolution (see also fig. [T|). Including source evolution, we find H = 101, 69, and 136, respectively. In deriving these 
estimates we have approximated the neutrino spectra as tracing the proton spectra; for softer spectra the results are 
somewhat lower. Scaling the neutrino fluxes presented in fig. 1 of Ref. [27j with these values, we find that the diffuse 
neutrino fluxes corresponding to model (ii) are well above the existing AMANDA-II limits [2{| [3(| in the case of strong 
source evolution. For all models except the linear accelerator (iii) the expected event rates in IceCube (obtained by 
scaling the results presented in table 1 of Ref. (27|) are 3> 1 per year, even in the case of no source evolution. Hence 
IceCube should be able to put strong constraints on these models. 
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IV. SUMMARY AND DISCUSSION 

In this work we have considered the relation between the diffuse neutrino flux and the neutrino flux from Cen A 
under the hypothesis that Cen A is a characteristic source of UHECR protons and UHECR-related neutrinos. This 
is motivated by recent results from PAO GlBl which suggest that Cen A may be a source of UHECRs. If it is also a 
source of neutrinos, as proposed in Refs. [26l. 127] . and if the environment in Cen A is representative for other sources, 
we argue that the diffuse neutrino flux may also constrain the source model. The diffuse neutrino flux can easily be 
estimated by scaling a model Cen-A neutrino flux (see eq. ((HJ) if one assumes that Cen A is a 'typical' source, i.e. 
under the assumption of universal UHECR and neutrino injection spectra with a fixed relative strength. We stress 
that we make no assumptions on the intrinsic luminosity or distance of the sources. We have derived estimates for 
the corresponding scaling factor in section [ill without relying on a particular source model (see fig. [I]). The scaling 
factor depends mildly on both the assumed proton and neutrino energy spectra, but very strongly on the assumed 
model of source evolution. This suggests that constraints related to the diffuse neutrino flux may be especially useful 
in constraining the evolution of UHECR sources. 

Regardless of the source model, we find that the estimated neutrino event rate in IceCube due to the diffuse neutrino 
flux is expected to be at least two orders of magnitude larger than the event rate of neutrinos from Cen A, as may 
be seen from eqs. © and (fTU]) . When sources follow strong AGN evolution as parameterized in Ref. [Hj], the rate is 
even three orders of magnitude higher. Therefore we conclude that the detection of neutrinos with IceCube, without 
the detection of the associated diffuse neutrino flux, would imply that Cen A is an exceptionally efficient neutrino 
source. Here we assume that the neutrino flux extends to energies above 10 6 GeV, so that neutrino detection is not 
limited by background rejection. 

We have applied our results to models recently proposed in the literature in section [TTT1 We find that the diffuse 
neutrino flux associated with the model adopted in Ref. [26[ is well above the (preliminary) AMANDA-II UHE limit 
(30| and the PAO limit [3l[ when the sources follow strong evolution. This implies that either the model (when applied 
to strongly evolving sources) overpredicts the neutrino flux, Cen A is intrinsically exceptional, or that considerably 
more than two of the UHECR events observed by PAO are produced by Cen A. Similarly, we find that the diffuse 
flux associated with the most optimistic models considered in Ref. 27] is also above AMANDA-II limits [1^, [3(| for 
stron g so urce evolution. The expected event rate in IceCube is much larger than 1 yr _1 for all models considered in 
Ref. |23| except for the linear accelerator (both in the case of no source evolution and of strong source evolution). 
IceCube should thus be able to put strong constraints on these models. 

Several comments are in order. First of all, it is presently unknown to which extent UHECR source are similar 
in nature, and whether or not Cen A is a typical cosmic-ray source (recall that we use 'typical' only as a statement 
on the ratio of outputted neutrino flux to cosmic-ray flux; it bears no meaning to the intrinsic luminosity or the 
distance). We have found no a priori reason to suppose that Cen A is an atypical cosmic-ray source. Cen A is quite 
representative of the local population of radio- loud AGNs 0] , which, although constituting a subdominant fraction 
of 15 — 20% of all AGNs [42|], are commonly (though not uniquely) considered as the dominant cosmic-ray sources 
because of their powerful jets. Notwithstanding this plausibility argument, the question whether or not Cen A is 
indeed a typical UHECR source can only be resolved by further observations. Our results address this question by 
placing constraints on the simplest scenario, namely that Cen A is a typical source, representative of a universal class 
of cosmic-ray sources. In this context we reiterate our conclusion that the detection of neutrinos from Cen A, without 
the detection of the diffuse neutrino flux, would rule out Cen A as a typical cosmic-ray source on the ground of its 
exceptionally high neutrino production efficiency. 

A second comment regards the number of UHECR events attributed to Cen A. In our numerical estimates we have 
followed the PAO analysis, which associates two events with Cen A as the angles between their arrival directions 
and the nucleus of Cen A are smaller than 3.1°. It has been pointed out in Ref. Q that at least four events can be 
associated with Cen A if one takes account of its extended structure (the radio lobes subtend on the sky) . On 
the other hand, it is also possible that the observed UHECR events are not from Cen A but rather from sources in 
the more distant Centaurus supercluster, as suggested in Ref. In both cases our estimates on the neutrino boost 
factor would be reduced by a factor ~2. This would not qualitatively affect the conclusions for the models studied in 
this work. Future data from PAO will increase the rather limited statistics and will very likely shed more light on the 
UHECR production rate of Cen A. 

Thirdly, the relation between the produced neutrino and UHECR fluxes may be complicated when neutrinos and 
protons are emitted within cones of different opening angle (as in the model discussed in Ref. [43]), or when the 
source luminosity varies in time. These effects will average out for diffuse fluxes but may affect the emission from 
Cen A. Any collimation in the direction of the jet will reduce the visibility on Earth as Cen A is viewed off-axis. 
Because the opening angle of the neutrino emission cone is expected to be smaller than that of the proton emission 
cone, this further decreases the expected neutrino flux from Cen A compared to the diffuse neutrino flux. Radio and 
X-ray measurements of Cen A show variability on a time scale shorter than a year Q. The arrival times of protons 
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and neutrinos produced during a flaring phase will however not be correlated as the proton path length is increased 
by its motion in the galactic magnetic field. Hence a strong flare may in principle lead to an increased neutrino flux 
without an increase in the observed UHE proton flux from Cen A. 

Finally, we note that the diffuse gamma-ray flux produced by sources similar to Cen A may also constrain the 
source model, in a fashion very similar to the diffuse neutrino flux that was considered here. Gamma-ray emission by 
Cen A was considered in Refs. [23, E3 but these authors have not considered the associated diffuse gamma-ray flux. 
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APPENDIX A: PROTON ENERGY LOSS APPROXIMATION 

Proton energy loss during propagation from a source at redshift z can be described by the differential equation (see, 
e.g., Ref. [H) 

IdE 1 (l + z)(3 ((l + z)E) 

E dz 1 + 2 H(z) ' 1 ' 

where the first term accounts for redshift energy loss and the second for particle interactions. The function @q gives 
the inverse energy-loss time (at present epoch) for interactions between protons and CMB photons. 

In this work we determine the energy at the source Eq as a function of observed energy E and redshift z by solving 
cq. (| Al|) numerically. To do this, we split /3q = /3J + /3q £ into two parts corresponding to energy loss due to pion 
photoproduction and electron-positron pair production, which we approximate by: 

5 

log/3J = Y.anX*- 1 ; (A2) 

n=l 
5 

log/3 ee = ^2b n X n -\ (A3) 



where X — log(£'/leV), and (3$ is expressed in units of yr 1 . For 10 10 5 GeV < E < 10 12 GeV, a = 
(—1.2 x 10 5 , 2.3 x 10 4 , —1.7 x 10 3 , 52, —0.62); for higher energies a\ = — 7.6 is the only non-zero coefficient; for 

lower energies /3J = 0. For 10 8 5 GeV < E < 10 12 GeV, b = (-1.3 x 10 4 , 2.6 x 10 3 , -2.0 x 10 2 , 6.6, -0.082); at 
lower or higher energies (3$ e = 0. 



APPENDIX B: NEUTRINO DETECTION WITH ICECUBE 



In this section we estimate the effective area for the detection of muon-neutrinos with IceCube. The word 'neutrino' 
refers to muon-neutrino throughout this section. 



1. Neutrino interaction and muon track length 

The expected number of neutrino events in IceCube for a neutrino source with differential flux (j^iEy) at an angle 
9 with respect to the nadir (i.e., 9 — points towards the North Pole) can be written as follows: 

N = T J dE v <p*{E v )V v {E v ,6) J dEp n(E u , E^T^E^A^.^E^) , (Bl) 

where T is the observation time; E v the neutrino energy; E^ the muon energy; V v the probability that a neutrino 
reaches the vicinity of the detector; the probability that a muon is created that reaches the detector with sufficient 
energy for detection; A^ )e g the detector's muon effective area (which is close to the geometrical surface for high-energy 
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muons); and n(E v , E^) is the muon energy distribution resulting from the interaction of a neutrino with energy E v . 
In the following we evaluate these quantities and make a number of simplifying assumptions. 

Neutrinos interact with Earth nuclei through both charged- and neutral-current interactions. Interactions of the 
former type lead to electron, muon, and tau production while interactions of the latter type degrade the neutrino 
energy. Here we assume that all interactions transform the neutrino and thus we neglect the 'regeneration' of lower- 
energy neutrinos by neutral-current interactions. The neutrino survival probability V v can then be expressed as 

V V {E V1 6) = cxp ^-N A a vN {E v ) jf p{r)dl , (B2) 

where N A = 6.2 x 10 -23 cm -3 water equivalent (w.e.) is Avogadro's constant; er„jv is the total neutrino-nucleus cross 
section (including charged-current and neutral-current interactions) ; p(r) is the density of the Earth as a function of 
the radial coordinate r = \Jl 2 + r' E — 2We cos 9; and L = Ite cos 9 is the propagation distance for a neutrino through 
the Earth with nadir angle 9 (te — 6.4 x 10 8 cm is the radius of the Earth). Since there is no experimental data on 
neutrino-nucleus interactions at the UHECR energy scale one has to rely on models to extrapolate the data from lower 
energies. In this work we use the neutrino-nucleus cross sections tabulated in Ref. [451 ]. We adopt the Preliminary 
Reference Earth Model (see Ref. (45|) for the density profile of the Earth p(r). The muon detection probability 
is: 

V^Ep) = 1 - cxp (-N A a^(E v )R^{E^)) , (B3) 

where er^jy is the charged-current cross section, and is the muon range within which the muon energy degrades to 
a minimum energy E^ m . This range can be approximated with 



, , 1 / a + bE u \ , s 

Ra(Eu) = T n — A- , B4 

where a = 2.0 x 10 _3 GcVcm _1 (w.e.) accounts for ionization and b = 3.9 x 10 _6 cm _1 (w.e.) for radiation losses 
[IB] ]. In this work we adopt £™ m = 10 2 GeV. For simplicity we assume that a neutrino interaction leads to a single 
muon with energy E^ = ycc(E v ) E v , i.e. 

n(E u ,E fi )=5(E fl -y C c(E u )E v ), (B5) 

where the charged-current inelasticity ycc is tabulated in Ref. [45| . Lastly, we use the IceCube muon effective area 
given in Ref. [39| ■ This is the only quantity in our estimates that accounts for detector efficiency. 



2. Effective area for downgoing neutrinos 



Since downgoing neutrinos reach the detector virtually unhindered, we may approximate V v ~ 1 and ~ 
Na<J^(E u ) mm (R^, Rd), where Rd — 2 x 10 5 cm denotes the detector depth. At the energies where downgoing 
muons are detectable the interaction length is determined by the detector depth (i.e., Rd < Rfj,)- With these simpli- 
fications, we express the expected number of neutrino interactions for a single point source above IceCube as 

N^ dn =T j ' dE v <t>f{E v )A%s{Eu) , (B6) 

where T is the observation time, is the neutrino flux from the source and A dn cS denotes the neutrino effective area 
for downgoing neutrinos: 

A% S {E„) = TRdN A a^{E u )A^ME^ ■ (B7) 

For simplicity we take the muon effective area for downgoing neutrinos equal to the geometrical area, A^^g = 10 10 
cm 2 , so that the corresponding neutrino effective area is angle-independent. In this case we may express the number 
of events due to the diffuse flux of downgoing neutrinos as follows: 

^diff.dn = TQi J ^ ^^ Ev)A An^ Ev) t (Bg) 

where fii is the detector's opening angle and 4>f, lS is the diffuse neutrino flux. We take the maximum viewing angle 
for upgoing (downgoing) events equal to 5° below (above) the horizon, so that = 5.7. 

We note that the detection of downgoing neutrinos is challenging and requires special analysis techniques. Conse- 
quently our estimates for the detection rate of downgoing neutrinos may be too optimistic. 
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FIG. 3: Estimated neutrino effective area for the detection of upgoing and downgoing muon neutrinos with IceCube. For 
upgoing neutrinos, we have chosen iJ™ m = 10 2 GeV. For downgoing neutrinos, we adopt the geometrical muon effective area 

A^e = 10 10 cm 2 . 



3. Effective area for diffuse upgoing neutrinos 

For upgoing neutrinos the event rate in angle-dependent because the neutrino survival probability depends on the 
incident angle. We compute the number of neutrino events for the diffuse flux of upgoing neutrinos by integrating eq. 
(|Bip over the angle 9. We express the result as follows: 

N diS, up = rfii J ^ 0diff ( ^ )A uP a(£ g _ (B9) 

Here the average effective area for diffuse up-going neutrinos, using the simplifying assumptions described above, is 
given by: 

AZ S (E„) = S(EME^ME,) , (BIO) 
where we adopt values for A fl _ c ff(E^) from Rcf. [3^ |. and the shadowing factor S is: 

S(E V ) = - / d6am6Vv(E vt 6). (Bll) 

1 - cos max J 

Here # max = 85°, i.e. 5° below the horizon. 

In figure [3] we show our estimates of the IceCube effective area for downgoing and upgoing neutrinos, given in eqs. 
(|B7|) and (|B10j) , respectively. The initial increase of the effective area for upgoing neutrinos with energy is due to the 
increased muon path length, whereas the subsequent decrease is due to the fact that the Earth becomes opaque to 
neutrinos. These effects play no role for downgoing neutrinos; here the energy dependence of the effective area follows 
the energy dependence of the neutrino interaction cross section. We stress that the effective area for downgoing 
neutrinos corresponds to a muon effective area A^^g = 10 10 cm 2 . Given the difficulty in detecting downgoing 
neutrinos, this is very optimistic and may be considered as an upper limit. 

We have verified that our estimates on the neutrino effective area for the upgoing diffuse flux agree with the 
estimate presented in Ref. (4(| within a factor two (for energies below 10 7 GeV) to three (at 10 s GeV). Furthermore, 
the corresponding event rates for a fiducial <j) v cx E~ 2 source spectrum are within a factor 2 of results presented in 
Ref. [H. 
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